

Here we present new calculations of optimal parameters for Uda-Yagi ** 
plane antenna with elements parallel to ideally conducting Earth. 
Antenna consists of 6 cylindrical elements diameter 2a each, namely 
reflector, emitter and four directors with lengths slightly varying to give 
the desirable radiation pattern . 
Calculations are aimed at optimizing of parameters 
( maximizing forward/back intensity ratio called "shield coefficient 
and/or direct wave intensity) 

We're calculating current distributions in the reflector, emitter and 
all the directors which provide zero electric field on the side 
surfaces of them and turn into zeros at the edges. In other words 
we find a numerical solution of boundary- value problem for 
Maxwell equations, Zommerfeld's radiation condition included, 
for the Uda-Yagi antenna. 

For all the numerical results hold current symmetries discussed 
below and reciprocity theorem. 

Antenna parameters, i.e. lengths of elements and mutual distances 
are adjusted during iterations process , first iteration being 
the results of MMANA program implementation to the antenna . 
Our calculations showed that the MMANA program applied to 
the Uda-Yagi antenna provides a system with forward intensity 
amplification of 25 times compared to single cylindrical antenna , 
both at height of 4 meters above ideal Earth, with electromotive 




force (later EMF ) of external feed in the middle of emitter the same. 
Shield coefficient (i. e. forward/back amplitude ratio) amounts to 
20 decibels for Uda-Yagi antenna with MMANA parameters and 
total intensity ratio for Uda-Yagi to single antenna approximately to 3 . 
Our iterative adjustment of lengths gave direct wave enhancement of 40 
instead of 25, shield coefficient of approximately 49 decibels instead 
of 20 and total intensity ratio of slightly more than 4. 
We split the reflector, emitter and directors into sets of infinitesimal 
current elements (dipoles) and calculate total longitudinal electric 
fields on the side surfaces of all such elements through spacings of 
elements with unknown infinitesimal currents. The number of 
infinitesimal currents to calculate is obviously equal to the number 
of elements with zero longitudinal field on side surface 
So we arrive at the system of linear equations for unknown complex 
amplitudes of infinitesimal currents, having single solution , best 
represented in tensor form. 

Tensor form of system of linear equations enables us to implement 
for the solution a procedure quite analogous to Gauss-Jordan 
elimination or Gaussian elimination with back substitution. 
We formulated our linear system so as the latter appeared 
transparent and slightly easier and without any noticeable rounding 
errors in matrix calculations with 1 6-digits computer accuracy. 
Linear system consists of 6 tensor equations, one for reflector, 
emitter and four directors, with 6 vectors of infinitesimal complex 





amplitudes of current elements (dipoles) to be determined and 
6*5=30 matrices of coefficients calculating longitudinal electric 
fields or, more precisely, the distributions of induced local EMF. 
Currents produced by these electromotive forces produce fields 
zeroing on side boundaries of each element of Uda- Yagi antenna. 
Multiplication of coefficients in Gauss elimination is plainly 
replaced with matrix multiplication and division of coefficients - 
with calculating a reverse matrix. All matrices to be reversed 
appear to be perfectly conditioned ones so we did not need the 
regularization procedure. Suffice it to say that, Crout method 
with or without pivoting gave no difference. The solution thus 
obtained fitted the linear system with accuracy of 16 digits. 
We now replace infinitesimal current elements with discrete ones 
with lengths determined by their number on each antenna 
element (chosen equal to some constant N for all elements i.e. 
reflector, emitter and four directors) 

The i-th current element complex amplitude is equal to the Green's 
function (to be determined ) of antenna element multiplied by 
local EMF with summation over all current elements (summation 
over index k ) from 1 to N. 

h = g, k e k (0) 

EMF in the point k on antenna element "m" created by all current 
elements of antenna element "n" ("m" and "n" denote any pair 





of antenna elements) equals to 

e k =f,Ij (1) 

EMF matrix /,. here represents EMF (i.e. field strength multiplied 

by discrete current element length ) in the point k of antenna 
element "m" created by the current element in the point j of antenna 
element "n". To get the EMF in the point k of antenna element " m" 
by all the current elements on antenna element " n " the summation 
over index j is required. 

From (0) and (1) taking an account of the associativity of matrix 
multiplication we see (" m " and " n " are notations of any pair 
of antenna elements ) 

I m >={g,J kj )I" j (2) 

Value in round brackets represent current influence matrix of 
antenna element " n " imposed on antenna element " m " 

Let us use the notation for antenna elements : r - reflector , 

e- emitter (i.e. center- fed antenna element) , f- first director , 
s- second director, t- third director, 1- (last) - fourth director 
and two-letter names for current influence matrices-first letter 

denoting antenna element " m " in (2) and second letter 

denoting antenna element " n " in it. 

Here we readily recall that in the matrix product ( g ik f kj ) first 

matrix is the antenna element " m " Green's function depending 
of element " m " only. 





and second matrix is the mutual EMF matrix depending of mutual 

spacings and properties of both. 

The reciprocity theorem states that when antenna elements m and n 

exchange, EMF matrix can be obtained by transposition. 

Albeit this holds , we have no need to use it. 

In terms of the notations taken, for example 

&,=(*•/*) (3) 

here lt t . is influence matrix representing part of current in i-th point 

of fourth (last) director induced by current in j-th point of third director 
First matrix in round brackets is Green's function of last director and 
second matrix is the mutual EMF matrix of both. In view of reciprocity 
theorem, following same notations, tl tJ = (g ik q kJ ) here q kj = f jk , 

because Green's function g ik obeys the reciprocity theorem, too 

Let us find an expression for EMF matrix. Suppose we neglect any 
azimuthal currents, so vector potential has only one component 
parallel to antenna element. 

For the vector potential of current element, in obvious notations, 
dA = I*dl /(c*r) , and for electric field vector 

dE = -—(graddivA + k 2 A), [2] (4) 

ik 

Let's place the center of polar coordinates in the center of current 

element with polar axis parallel to the element of antenna, then. 

(u is the polar angle) 





divA = cosu(-^- + — ) — — — -I* hi 
r r c 

After simple calculations we get (h2 is the length of radiating current 




element) E r = 1 cos v — — - (- + —^ 



I* hi exp(ikr) / 1 i 
V kr 1 



i 1 exp(/A:r) 

kr r cr 



E v =smv(^ Y + --ik)^^I*hl 



Longitudinal field E = E r cosu - E v sin v , that is 

E, = h * h2 ex P(^) (( 1 + -i-) * (3 cos 2 u - 1) + ik sin 2 u) 
c r r kr 

The EMF in any current element of any other element of antenna if 

the product of field by length of that other current element, hi. 

e=f ;J = Ij * h2 ew(ikr) ((- + ^-)*(3cos 2 v-l) + iksm 2 v)* hi (5) 

c r r kr 

with r and v determined by indexes i and j . 
Complex amplitudes L are values to be determined 



Next we come to the task of calculating the Green's function 
for antenna element parallel to ideal flat Earth at the height H, 
given the Green's function of that antenna element in the free space. 
Boundary conditions hold if image antenna at same 
" height " H below Earth and of opposite current direction added [2] 
Let us seek for the Green's function in the form 

Sik = g°% + x ik ( 6 ) 

with Green's function to be found in the LHS and given Green's 
function for free space and appendix x ik ( not necessarily small) 




caused by field of image antenna. 

According to what has been said we get 

% = Ry (rg°fi - Xjk ) where R tj = g°afy 

Matrix in round brackets is the Green's function of image 

antenna and f« - EMF matrix of image and the original antenna. 

From the last equations we derive ( 5 tJ is Kronecker ' s symbol ) 

(Ry+S^xjt =-R ilg ° lk S y x jk =x lk =-(R lJ +S lj y 1 (R jl g° lk ) 

The matrix is round brackets in RHS of second equation is 
the reverse of that in round brackets in LHS of first equation. 

Here we come to the closed linear system of tensor equations with 6 
N-dimensional vectors of current complex amplitudes to be determined 

and 30 N * N matrices of influence as coefficients. First letter in the 
matrix name denotes the antenna element in which current distribution 
is calculated and the second letter denotes the antenna element the electric 

field of which produces this part of current . Each matrix is the Green's 
function matrix of first antenna element multiplied by the EMF matrix 
of first and second elements calculated above (3) 

Letters in matrix names mean r - reflector, e- center-fed emitter,f- first 
director, s- second director, t- third director, 1 - the last (fourth) director 

// ={lt) tJ I) +{ls) iJ I) +(1/1. If +(lr) tJ I r j +{le) tJ i; (7) 

• l\ ={tl) lJ I , J +{ts) ij I)+{tf) ij I^{tr) ij r j+ {te) ij i: (8) 





• i: =(sl) i jIj+(st) ij I< + ( S f) ij lf + ( S r) ij r j+ ( S e) iJ i; (9) 

• i{=(fi\ji l j +(/t) tJ i t j +(/4/; +(|/>),/; +(fe\ji; o°) 

• ^=(rl) 1J I-+(rt) iJ Ii+(r S ) ij I° + (rf) ij l{+(re) ij i: (11) 

. /; =/°, +( e o„/i +(4^J +M^y +M// +(«4 A' (12) 
The additional current 7°, in equation 12 represent the KNOWN current 
produced by external EMF exerted over the emitter in the center of it. 

The ideal Earth's influence is twofold: through the Green's function modified 
as above and additional EMF analogous to that calculated in (7) caused 
by image antenna. Both corrections are to be subtracted from 
original expressions cause of opposite current in the image antenna 

Here we are at last at the crucial point of calculating the Green's function 
of a single conducting cylinder of length 2L and diameter 2a in free space. 
Notations taken are shown below: cylindrical coordinate system with Z axis 
being the axis of cylinder and zero (center) in the middle of cylinder axis. 

7? = -yj(z-z'} 2 +a 2 +r 2 - 2ar cos(q> - q>') - the distance between two points, one 

of which with coordinates z' , q>' is on the cylinder's side surface 








Let us consider an integral over the cylinder's side surface 



In L 



Qxp(ikR ) 



V(z)=\\ — — -S(cp',z') a dq>' dz' S = J/c j- surface current density 



-L 



R 



. V(z)- vector potential, f = exp(ikR)/R - Green's function of the 
Helmholtz equation which takes the form ( k= In IX - wave number ) 



I d( t df\ t d 2 f 1 d 2 f 
dz 2 + r 2 dcp 2 



r dr{ dr 



k 2 f = 



V vrj 



Separation of variables entail ( w and s are parameters ) 
f = R{r)Z{zyP(cp) ±^ + w 2 =0 -^ + s 2 =0 



Z d z 



®dq>' 



Radial function R(r) obeys the equation 



r 



.2\ 



w l -k l +■ 



R = 



J 



It is limited on the z axis ( r = ) and tends to zero at infinity. 




The Green's function of Helmholtz equation can be expanded into the series 




JkR 



- — = ^s s f s cos(s((p - (p')) with s = 1 , £j = s 2 = • • • = 2 Neumann coefficients 
R s 



In ikR 



1 a 1 ' 

chosen so as we should have f s = — [ cos(s(<p - <p')) d(<p-<p') 

ln R 



Let us neglect any azimuthal variations, in compliance to above . 
Thus s = , and the equation for R(r) becomes the modified Bessel equation 
of zero order and all Fourier coefficients except for f (r, z) denoted further 

as f(r, z) , become zeroes. Separation of variables thus give 



/ = A jcosw(z-z') 



I W ir " 



k l r) K n Ww-k z a 



K [Jw 2 -k 2 r)l o y w 2 -k 2 



dw (13-0) 



where the modified Bessel solutions, i.e. Infeld function I and 



Mc Donald function K are chosen so as f(r, z) be continuous at 
r = a. A - a coefficient to be determined. Upper line in figure 
brackets should be taken when r < a lower line - when r> a. 
Comparing f(z- z' , r ) to the vector potential in the point z = z = r = 
created by dipoles uniformly distributed over the cylinder's basic circle 



(Appendix 2 ) we get A = — , so finally 

n 



1 °f 
f = — \ cos w[z-z) 
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IqU/h 



2 1,2 



k l r K a y W z -k 2 a 



K [Jw 2 -k 2 r) I [Jw 2 -k 2 a 



dw 



(13) 



This expression is exactly the averaged fundamental solution exp(ikR) / R 




in the terms of reverse Fourier transform of the value in figure brackets. 
Along with the distance R, it is an even function of the difference 

of z coordinates. This enable us with the possibility to find the 
current distribution U(z) and hence the vector potential V(z) in the 
form of Fourier series, which we treat as " Least Mean Square Deviation" 

series , taking no account of Dirichlet theorem (3) 

When w < k we have vw -k 2 = -i-Jk 2 -w 2 \yjk 2 -w 2 > Oj , 

in view of Zommerfeld's radiation condition . Taking an account of 

I (-ix) = J (x) K (-ix) = i— H„(jc) if w < k , we get when (r<a) 




I {jw 2 -k 2 r) K Q {jw 2 -k 2 a)=i-J {jk 2 -w 2 r) U l (4k 



2 -w 2 a 



and when 

(r>a): I [Jw 2 -k 2 a) K [Jw 2 -k 2 r)=i-J []k 2 -w 2 a) U\(yjk 2 -w 2 r) 

The above integral representation for vector potential V(z) 



Jf— ^ -S(<p',z') a dq>' dz' -V{z) S = J/c j- surface current density 



-I 



with the not depended of azimuthal angle U(z) = S = I(z) / c 
in left-hand side with I(z) = 2^aj(z) becomes 

L 

\f(z-z',a)U(z')dz' = V(z) (14) 

-L 

with integrating along the cylinder axis length. 

We will further omit here cylinder radius a as constant and consider this 
as an integral equation for current distribution U(z') . 





Let us write the vector potential as follows 



V(z) = C°V°(z) + C cos kz + C'sin kz 



(15) 



The first term in the right hand side represents vector potential in the 
point with coordinate z caused by the current induced by the EMF 
applied in the point with coordinate z equal to unity produced 

by some external field shaped as Dirac delta function , 

E(z) = 5(z - z ) . It is being calculated easily, on account of three 

simple relationships to follow. 

First expression represents the zero balance between the external 
delta-shaped field and the field of current induced by external EMF 

equal to unity (integral of the delta- function) . Second expression 
relates the induced field to the vector potential to be defined. Integrating 
in the vicinity of point z we get the third expression 



E z +Et = E = -{ 



i ,d V 



z z 



k dz 



+ k z V" ) 



dV ( 



dz 



_dr 

z = z,, + dz 



z = z -0 



= ki 



after which we find for the vector potential easily 



V°(z) = -smk\ z-zA C°=l 

2 ' ° ' 



(16) 



The independent solution of this, V (z) = 0.5 * exp(zl| z- z |) leads to 
the very same current distribution, as we shall see shortly afterwards. 



It is essential for the technique to follow, that along with classical Fourier 
series in sin(nx) ,cos(nx), "half-integer" Fourier harmonics 




sin((n+0.5)x) cos((n+0.5)x) obey the orthogonality condition on 
the same interval .We draw this conclusion on the basis of formulas 
quite analogous to these describing the Fourier expansion outlined in [3]. 
Hence the half-integer Fourier harmonics can also be used for an expansion 
analogous to Fourier series with no zero term and Neumann factor equal to 2. 
The kernel of integral equation (14) is even function depending on the module 
of its argument , so we can divide current distribution and vector potential 
into even and odd parts and get the solution of (14) for both independently. 
For even part we use the ordinary Fourier expansion and for odd part 
- an expansion in half-integer Fourier-like series. Period of f(z - z')) 
equals 4L, periods of U(z) h V(z) equal 2L. 
. Let us first consider the even part of solution and write the Fourier 

expansions of the vector potential in the RHS , kernel of (14) and 
current distribution in the form 




2L 



F v = \f(z) cos vgzdz (17) 



/= ^s v F v cos vg{z-z') g = ^- (18) 



CO CO 

V(z) = Y,£ ll V, 1 cosv -z = ^s n V n cosngz (19) 

77=0 L n=0 

U(z) = ^s q U q cos( qgz) e =l £ 2 =£ 4 =... = 2 (20) 
9 

Certain asymmetry in (18) compared to (19) and (20) is made for 

some future convenience. Values of n and q in (19) and (20) are 

even ones only. This is because of that U(z) and V(z) periods are 

a half of f(z) one. 





On account of formula (13) and its reliance to Fourier transform of 
value in figure brackets there, taking an account of (17) we obtain . 

K =I {j(v 2 g 2 -k 2 )a)K ^(v 2 g 2 -k 2 )a)+AF v AF V =-] f(z)cos(vgz)d z (21) 



21 



Modified Bessel functions were calculated with 1 6 digits computer accuracy 

in accordance with algorithms in [1] . 

In calculation of the last integral in (21) we can make use of an approximation 

e ikz 
f(z) = — with simple corrections albeit we now bear in mind some 

z 



much more sophisticated procedure. 
[ cosvgzdz cosvgz =— (e' vgz +e~' vgz ) 



2L Z 



°° „iz(k-¥Vg) °0 fg 

[—7- ^du= f —d% 

i L z{k + vg) 2k l vji B, 

As a result 

i oo it 

kF v =-{E{2kL+vn) + {2kL-vn)) E(x)=Ci(x)+iSi(x)=-\ — dt (x>6) 

X 

E(x) = E* (-x) = Ci(- x)-i Si{- x) (x<0) 

1 L 
The Fourier theorem gives V n = — f v(z) cos ngz dz (22) 

and after substitution of (14) 

1 L 
V n = — \\f{z-z')cosngz U{z')dzdz' (23) 

After substitution of the U(z) in the form of Fourier series (20) 




we obtain the system of linear equations with coefficients D nq : 

V„=J J D nq U q ;D nq = ^-\\f{z-z')cos(ngz) cos{qgz') dz dz' (24) 

The integral of D nq can be calculated easily by the substitution of series 
(18) for f(z-z') with subsequent term-by-term integration. 

nq ~ tfinq ' 2 q / ,nm &am "m 

n m 

r , \ l,T, " TI ffi 

a nm = g cos (ngz) cos(mgz)dz = {-\) 2 — (26) 

* n -m 

Here n and q are even values, m and p - odd ones. 

Details are quite obvious: we have 

g « 
/ = — 2^ £ v F v ( cos ( v S z ) cos(vgz') + sin( vqz) sin( vgz')) 






(2 


5) where 


n+tn+l 

2 




m 


n 2 


-m 2 



and double integral turns into multiplication of two simple ones ,with 
results proportional to a nm 

Obviously, all integrals of sin (ngz) * cos (qgz) and vice versa are equal to zero 

and integrals of cos (ngz) * cos (qgz) are nonzero if n = q only , and so we 

come to the formula (25) . 

For the even part of V(z) we have 

V(z) = C°V°(z) + C cos kz 

The Fourier transform of second is quite simple, and so we have 

V„=CY° + CV: K = ^ V'=AkL sin kL- ^ 



kL n 2 (2kL^ 

V n J 



2 

n 




On account of (16) 




V°(z) = — smk\ z-z n 
2 ' 



1 r 



V\ = — j" V°(z) cos ngzdz (27) 



After substitution of first of (27) into the second and integration we get 



K ( z o) = 



4kL y 



.2 



K 



cos((^L + ?7^:)x ) co%i^kL-r\7i)x^) 



2kL/Y _ n 2 (k + ng)2L (k-ng)2L 



cos 



(ngz ) 



+ 



+ ■ 



in (ngz ) 



sin 



sin(kL+rj7z)x sin((kL-rj7z)x ) 
(k + ng)2L (k-ng)2L 



+ 



+ ■ 



cos 



{ngz ) 



AkL, 



K 



2kL, 



n 



cos{i y kL + rj7z)x l ) cos^kL-rj^x^ 
(k + ng)2L (k-ng)2L 



sin 



(ngz ) 



sin^ikL + rj^x^ sin [kL-r]7r)x 1 



(k + ng)2L (k-ng)2L 



And finally, for even part of vector potential 

V n °(z ) =0.5*(Vj (z ) + V; (-z )) 

The second solution of (16) , namely, V 00 (z) = 0.5 *Qxp(ik\ z-z |) 

differs from the first one by an item of 

, V 01 (z) = 0.5 * (cos(k * (z - z )) + cos(£ *(z + z ))) = cos(kz) * cos(kz ) 

Hence the implementation of it differs constant C in vector potential 

V(z) = C°V°{z) + C coskz 



only, and not the current distribution, because k and z are constants 





Here x =l-^- x, = 1 + ^- r\= n A 

L 2 ' /2 

In view of linearity of Maxwell's equations and Fourier transform we 

may omit factor of i / 2 here and multiply the Green's function by it later 

The parameter z in the Fourier series equals that in Green's function G (z, z ) 

Indexes of g°,k (6) denote arguments of G (z, z ) with arbitrary order 

in view of reciprocity theorem. Boundary conditions give 

^A~¥u n =0 or U =-2^(-l)" A U n (28) 

n 

Here and later " shaded " sum means no zero term .For 2L periodic 
U(z) expansion in " half -integer " Fourier-like series turn into 
expansion into odd harmonics 
U{z) = ^s m U m cos mgz m=l,3,5.... (29) 

in 

quite analogous to the way the ordinary Fourier series is replaced by 
an expansion in even harmonics n or q U(z) = ^s„ U n cos( ngz) (20) 

n 

and coefficients U m ,U n obey a simple relationships 

1 L 
Analogous to U n = — \u(z) cos ngz dz (30) 



1 L 
for odd harmonics we have U m = — \u(z) cos mgzdz . (31) 

1L -L 

Substituting U(z) = ^s n U n cos( ngz) into (31) we get 



U =—Ya s U , (32) 

m / j n in n n ' V ' 




and substituting (29) into (30) we find U n =—Ya nM s m U m (33) 

n 




These relations along with some symmetrical ones to follow we need to 
establish the plentitude of Fourier- like series in odd harmonics. 
For odd parts of vector potential V(z) and current U(z) we write 



expansions in odd harmonics m and p 

v ( z ) = X s <» V <» sin m S z u ( z ) = H £ p U p sin (Pgz) 



(34) 



and, following the logic outlined in the calculations for even part we have 
a system of linear equations V m =^E mp U p (35) where 



o 

E „, P = F J mp + — £ Pnm Pnp S n F n (36) 

K n 



P nm = g J sin ngz sin mgz dz = (- 1) 



n+m+l 
2 



2 2 

n -m 



(37) 



with vector potential v(z) = C° V° (z) + C sin kz 



and Fourier coefficients 



V_ = C° V.! + C VI, VL =ArkL cos kL — ^ 



m + l 
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2kL 

V K ) 



m 




V n (z a ) = 



4kL, 



/ K 

(2kL/Y _ 

V /K> 



m 



cos((kL + /u7z)x ) cos((kL-/u7t)x ) 
(k + mg)2L (k-mg)2L 



sin 



(mgz ) 




cos 



(«gz ) 



sm([kL + ju7r)x ) sin((kL- /uk)x ) 
(k + mg)2L (k-mg)2L 



+ 



sin (mgz ) 



4kL, 



K 



2kL, 



K 



m 



cos((kL + idn)x x ) cos((^L-/i^)xj) 
(k + mg)2L (k-mg)2L 



+ 



+ 



cos 



(mgz ) 



sin((^L + y u^-)x 1 ) sin((AZ,- i u^')x 1 ) 
(k + mg)2L (k-mg)2L 



and, as a result, V m ° (z ) = 0.5*(F m (z )- F m (-z )) 

The second solution of (16) , namely , V 00 (z) = 0.5 * exp(ik\ z-z |) 

differs from the first one by an item of 

V 01 (z) = 0.5 * (cos(k * (z -z )) - cos(/c * (z + z ))) = sin(fe) * sin(fe ) 
Hence the implementation of it differs the constant C in vector 
potential V(z) = C°V°(z) + C sin(fe) 
only, and not the current distribution, because k and z are constants . 



here x n = 1 ■ 



L 



x, = 1 + - 
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M =- =(m-l)/2 + 0.5. 



Boundary condition (zero current at the ends) give ^(-l) 2 U m = 
Ordinary Fourier expansion of odd function U(z) is as follows 

1 L 
U(z) = ^s tt U tt sm(ngz)) Here U„= — j" U(z) sin ngz d z (38) 



and comparing this to the Fourier-like expansion in odd harmonics 





(second formula of 34 with m instead of p ) we find 

U m =~Y J s n p nm U n U n = -YPnm £ m U m (39) 

So ordinary Fourier series with the even harmonics due to choice 

TC 

g = — and Fourier like expansion in odd harmonics are 

equivalent., and by comparing them we get a sort of orthogonality relations 

\s Ys a a = \s Ts B B =5 (40) 

2 m /i n ^ nm np 1 m / > n r nm r np mp * ' 



\s Ys a a =\s Ys B B =5 (41) 

2 n ^^ m nm qm 2 « ^^ m * nm * qm nq \ I 

JL m JL m 



Putting (33) into (32) with replacement p->m and interchanging 
summations we get the first of ( 40) 

Similarly, putting second of (39) into the first we get the second of (40) 
First of (41) we get after putting (32) into(33) with replacement q->n 
and second of (39) put into the first gives the second of (41) the same way 

The characteristic of Uda- Yagi antenna is the so called " shield coefficient " 
i.e. the ratio of forward to backvard wave amplitudes. 
Our computations showed that for the accuracy of shield coeficient 
within 1 decibel the Fourier series of about 3000 harmonics should be 
put into calculations. 

Hence the implementation of the LU decomposition with the Crout method 
look dubious. . It occurs that iteration procedure leads much faster to solution 




A little transformation of (24) indicates the first iterations to employ 
For odd and even parts we have. 
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Orthogonality relations , used in the last equations of both, valid 
for infinite Fourier and odd Fourier-like series, are used to estimate 
the length of series required for the desired precision to attain. 
It appears that for the 8-digits precision we need finite Fourier series 
length of 5000 for F n and F m both, and, as final results showed, 

about 3000 -length Fourier series for " n " and " m " in 

V = IF U + Y d U V =2F U +Y e U 

n n n / , nq q m mm / , mp p 

1 P 

to attain the accuracy within 1 decibel for shield coefficient above. 
First iteration to take is that as if sums above are zero both, and we 
iterate these sums in both equations (series U n and U m to be defined) 

Constants C in expressions for vector potential provide the zero of both 

even and odd currents on the ends but now they prevent us from 

implementing the technique outlined directly . So we use the 

superposition principle and write current distribution as follows, 

for even part first : U„ = C° U° + CU' n 

and we get two linear systems with no undefined constants 
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First one represents current wave induced by EMF and second - reflected one 
Boundary conditions give us then 

IX(-i)^ 

C = -C° " 



I«.(-i)*u: 



where C° may be set equal to 1 with no loss of generality. 
For odd part of solution we have U m =C° U° m +CU' m 
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11+1 



S(-1)T-C/'. 

Solution procedure outlined converges to computer accuracy in 10-15 
iterations . Some other transformations of (24) are also possible and 

proven equivalent, for instance U n =C° U° + CU' n +U U n with zero 

even harmonic as a new variable, and " n " non-zero. Specifically, 
for even part we have 

V = C°V° + CV'o = D w U + X 'D 0q U q ;V n = 2F n U n + d n0 U + £ 'd nq U q 



here U a = C° U° a +CU' +U u" . Hence we obtain 



V n = C°V n ° + CV' n = 2F n (C°U°„ + CU' n + U U n ") + d n0 U + £ 'C% t/ c 



^ C <i B? £/', + V '[/„ <i B? C/" ? And come to equations for new variables 
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Zero harmonics, not used in these equations, are related as above: 
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+ CV'o (42) 



And boundary condition takes now the form 



cz ,/(-iyx +c/|z ;(-i) n/2 t/; 41 = - c *HS- i ) nl2u i ^) 



and these equations, after getting solutions of three linear system above, 
give us two constants C and U not defined yet. 



For some more ease in calculations we rewrite (42) as 
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(44) 



This alternative method obviously applies to the even part of current 
distribution only and gives, within computer accuracy of 1 5 digits, the 
very same results for both current distribution and vector potential 
(the solution of . of (43) and (44) is found immediately) . 
The calculation is completed by calculating the current distribution in 
the form of Fourier series for even part and Fourier - like series 
for odd, and multiplying the sum by a factor of 0.5 * i as noted above. 



External current I e in the emitter (12), created by EMF equal 1 is the 




Green's function with z = 0, and now we get a solution of (7)- (12) . 

This completes the task of finding current distributions and we 
go to field proceedings. The Green's functions found obey the 
reciprocity theorem with computer accuracy, currents calculated 
from (7)- (12) are zeros at the ends, and, with center fed emitter , 
current distributions appear to be symmetrical if antenna elements 
(emitter, reflector and directors) are center- aligned. 
The computer program based on this method completed all the 
calculations in less than an hour with 2.8 GHZ four-head processor. 
Six-head processor is expected to almost halve this time 
consumption as a little consideration shows immediately . 
For some practical Uda-Yagi antenna with elements chosen 
by MMANA program which gave shield coefficient of 20 decibels 
we succeeded in establishing shield coefficient of above 40 
decibels, by small corrections of the antenna elements lengths only. 

Appendix. . 

1 Let's calculate constant A in (13-0). For the sake of brevity 
assume k = r = . Infeld function equals 1 then and for z = z 

. CO 

/ = — = aA K (wa)dw We have an integral representation 

a 1 " 

u 




1 A 

K (x) = [ exp(-x * ch{t))dt Substituting into / = — = —[ K (x)dx 




and integrating 

-J CO 00 CO j CO j 

— =ff exp(-x*ch(t))dxdt= [ -^— = 2 f 7 ( y = exp(t) ) 

A J o J o { Ch(t) J (/+1) 

7r 2 

Last one equals — and 4) =— > accordingly. For symmetrical integration 

4 7T 

limits J =— . So (13) appears to be a reverse Fourier transform of even function 

n 

f(V-A) 




2 Let's give a Gauss elimination details of (7)-(12). 

It's common that equalities like A = A + B * C ( A, B, C- square matrices ) 

mean calculation of right-hand side and substitution the result into matrix A. 

E - identity matrix, i.e. Kronecker delta- symbol. 
First, we eliminate // in (8) substituting (7) instead : 
tr tJ = tr tJ + tl ik x lr kJ tf tJ = tf tj + tl ik + lf kj te tJ = te tj + tl ik x le k] 
ts tJ = ts tJ + tl ik x ls kJ tt tJ = tl ik x lt kj 
Now, transfer tt tj to the LHS and calculating the reverse of matrix there, 

a6 ij = {E jj - tt ij j ' , after multiplying matrices in the RHS by it 
tr tj = aO ik xtr kj te tj = a0 ik xte kJ t f tj = a0 ik xtf kJ ts ij = a0 ik xts kj 
we have instead of (8) 
// = (ts) IJ I) + (tf\. // +{tr) tJ I'j +(te) tJ I'j (8-a) 

Now, in the same way we eliminate // in (9) substituting (7) instead : 
sr t] = sr u + sl ik x lr kj se ij = se ij + sl ik x le kj sf tj = sf t] + sl ik x lf kj 




st tJ = st tJ + sl ik x lt kj ss ij = sl ik x ls kj 




and I\ , substituting (8-a) 



sr .. = sr .. + s t ik x tr kj se tJ = se tj + st ik xte kj 



Sfij = Sftj + st ik x tf kJ ss ij = ss iJ + st ik x ts 



kj 



then transfer ss tj to the LHS and after reversing the matrix there 
a6 t j = {E u -sSij) 1 after multiplying matrices in RHS by it 



sr tJ = a6 ik x sr kj se ij = a6 ik x se kj s f tJ = a6 ik x sf k 



kj 



we have instead of ( 9 ) 

/;=(,/), v //+(5r), v /; + (,4 v /; (9-a) 

As before, substituting into (10) equations (7) : 

fr,j = fr tJ + fl ik x lr kj fe tJ = fe ij + fl ik x le kj ft ij = ft iJ + fl ik x lt kj 

fs u = fs tJ + fl ik x ls kj ff tJ = fl ik x lf kJ , 

(8-a) : 

f r ij = f r u + fhk x tr kj fe u = fe u + ft a x *e kJ 



fstj = fs tJ + ft ik x ts kj ff tJ = ff tJ + ft ik x tj\ 



kj 



and (9-a) : 

f r ij = f r u + fs ik X sr kj f e ij = f e ij + f s ik x se kj ffij = ffij + f s ik x sftj 

and, after that, transferring ff tj to the RHS of (10) rewrite it as follows 

a9 u = ( E u ~ ffu Y l f r u = a0 ik x fhj fe tj = aO ik x fe kj we get 
l{={fr) i] I r ] +{fe) ij r i (10-a) 



Next, substituting into (11) equations (7) 




re tJ =re tJ +rl ik xle kj rf tJ =rf tJ +rl ik xlf kj r Sij =rs tj +rl ik xls 



kj 




rt = rt + rl . . x It, rr. ■ = rl , x /r, . 

'7 17 ik kj ij ik kj 



(8-a) 



™y =rs tJ +rt ik xts kj rf tJ =rf tJ + rt ik xtf k 



kj 



re u = re u + rt ik xte kJ rr tJ = rr tj + rt ik x tr kj 



(9-a) 



r fi i = r fi ,+ r s ik x sf k , re,. . = re, j +rs lk xse k . rr,. . = rr, + rs ,. k x sr k 



and (10-a) 



re, . = re, ,+rf lk xfe k , rr. . = rr,. t +rf ik x fr, 



kj 



rewrite (11) by LHS matrix reversion 

aG ij = ( E ij ~ rr u Y re ij = a0 ik x re kJ in the form 

i;={re\jl'j (11-a) 

and, finally, substituting into (12) equations (7) 



er tj = er tj + el ik x lr kj ef tJ = ef u + el ik x Z/ ly es iJ = es iJ + el ik x fa 
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ij = es tj + et ik x ts kj ef tJ = ef tJ + et ik x tf k 
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er . = er + et , x^r, . ee = ee +e^ -,xte, 



(9-a) 



e/y = e fij + es ik x */*; cr y = »ij + «*„ x sr /(7 ee, ; = ee,, + es. s x se 



/../ 



(10-a) 



er = er- ■ + e f- , x f r, ee = ee ■ + ef-, x /" e, 

'7 '7 J ik J kj ij U J ik J kj 




11-a : 

cc y =ee tJ +er ik xre kJ 

we obtain for the center- fed emitter 

(E-ee)..Ij =1° here E - unity matrix. 

Right-hand side represents current created by external EMF . Hence this 
equation for the total current in the emitter is solved immediately 

By means of back substitutions in equations (11-a) (10-a) (9-a) (8-a) and (7) 

we find all currents as follows. 

i;=re tj xl'j 

lf=fe tJ xI' J +fr tJ xI r J i; =se l] xl) + srij xl) + S f ij xlj 

I\ =te ij xl e . +tr ij xr. +t f iJ xl{+ts u xI s . 

I\ =le ij xl e j +lr tJ x/; +lf tJ xlf+htjxl'j +1^x1) 

Substitution this into original equations (7) - (12) showed that 
the solution obeys them with computer accuracy , appears symmetrical in case of 
center- fed emitter and elements center-aligned, as it should've been. 




Let us calculate the Fourier and Fourier- like transforms of vector potential (16) 
Rewriting the first one in the form we obtain. 

if. 1 z r 

V n = — I sin k(z-z Q ) cos ngz d z I sin k(z-z ) cos ngz d z;z-z Q =±x 

zr, -L 





V* n = — I sin Ax cos n g{x + z ) d x -\ I sin k(x) cos n g{x - z ) dx; 



On account of sin(kx) cos(ngx)=0.5*(sin(k+ng)x+sin(k-ng)x) 



and -sin(kx) sin(ngx)=0.5*(cos(k+ng)x-cos(k-ng)x) 



we have 
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Fourier -like expansion in odd harmonics appears to be 



If. . 1 z r . 
V m (z ) = — sin k(z — z ) sin mgz d z H sin k(z —z) sin m gz d z; z — z = ±x 

z -I 

1 ^ — z o 1 ^+ z o 

F*,«(z ) = — [sin Ax sin mg(x + z ) dx JsinA(x) sin mg(x-z )dx; 



and on account of sin(kx)cos(mgx)=0.5*(sin(k+mg)x+sin(k-mg)x) 



and -sin(kx)sin(mgx)=0.5*(cos(k+mg)x-cos(k-mg)x) 




we find 
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Uda-Yagi antenna image . r- reflector, e- center- fed 



emitter f , s, t, 1 - directors, from first to fourth (last). 



